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, Abstract. Massive stars in their late stages of evolution as Red Super- 

■ giants experience mass loss. The resulting winds show various degrees 

^0 ' of dynamical and chemical complexity and produce molecules and dust 

grains. This review summarises our knowledge of the molecular and 
^ • dust components of the wind of Red Supergiants, including VY CMa 

I \ and Betelgeuse. We discuss the synthesis of dust as a non equilibrium 

O . process in stellar winds, and present the current knowledge of the chem- 

istry involved in the formation of oxygen-rich dust such as silicates and 
^ . metal oxides. 
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1 Introduction 



The discovery of large masses of dust in high-redshift quasars and Damped-Lyman- 
Alpha Systems two decades ago (jPei et al. 19911 IBertoldi et al. 2003)) rekindled 
\ the study of cosmic dust formation in evolved circumstellar environments (ECEs) . 
' While evolved objects on the Asymptotic Giant Branch (AGE) are by far the 
" dominant dust makers in local galaxies and the Milky Way, other candidates con- 
tributing to the dust budget in the early universe include massive stars and their 
ultimate stage of evolution, their explosion as supernovae (SNe). The most abun- 
^ ■ dant dust-forming SNe are of Type II-P and have Red Supergiant stars (hereafter 
RSGs) as progenitors. The interaction of the SN explosion blast wave with the 
circumstellar wind of the RSG progenitor indirectly probes the nature of the RSG 
^ wind in term of chemical composition, dust content, and mass loss rate. Un- 
derstanding the way massive stars form dust thus becomes a crucial issue to the 
assessment of the galactic dust budget in the local and far Universe. 

RSGs correspond to the late stage of evolution of moderately massive stars 
(10 < M^, < 30 Mq). Theses stars show a wide spread in their mass loss rates, 
with values ranging from ^ 1 x 10~^ yr~^ for a Ori (jMauron fc Jos selin 20111) 
to ~ 3 X 10-"* Mq yr-i for VY CMa ( [Humphreys et al. 2005] ). The mechanisms 
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triggering mass loss are not well understood, and unlike the case of AGB stars, 
dust certainly plays a minor role to the overall mass loss process. RSG winds 
are probably accelerated in the stellar outer atmosphere by the combined ef- 
fect of convection, rotation, and chromospheric activity (jUitenbroek et al. 19981 
IJosselin fc Plez 20071 lOhnaka et al. 201 ip . resulting in an extremely heteroge- 
neous nebula where small to large convective cells, shocks, plumes and knots are 
observed (Kervella et al. 2009, 2011: [Chiavassa et al. 2010l) . Understanding the 
synthesis of solid grains in RSGs thus represents an interesting challenge because 
of the dynamics and harsh conditions encountered in the stellar outer atmosphere 
and wind. 



2 Molecules and dust in red supergiants 

A wealth of observational data exists on RSG winds in the (near-) infrared (IR). 
Molecules have been detected in the wind of several objects, including VY CMa, 
jj, Cep, and a Ori. The water molecule, H2O, was first detected in fi Cep and a 



Ori by Danielson et al. (1965) However, Tsuji (2000, 2006) first recognised that 
the observed near-IR absorption bands originated from a molecular layer called the 
MOLsphere, consisting of water, and to a minor extent, the cyano radical, CN. The 
MOLsphere is found to be located above but close to the stellar photosphere, at ~ 
2 R0. For the specific case of a Ori, a thin MOLsphere located at 1.3 R© made of 
silicon monoxide, SiO, water, and alumina (AI2O3) grains is necessary to reproduce 
the optical depth excess in the mid-IR (jPerrin et al. 2007)) . 

Major advances on the composition of RSG winds have been achieved with the 



submm Herschel space telescope. In VY CMa, Royer et al. (2010) showed that 
the wind harboured a very rich, oxygen chemistry dominated by H2O, where O- 
bearing species coexisted with carbon-bearing species such as HCN and GS. The 
presence of both oxygen- and carbon-bearing molecules may trace the breaking 
of CO molecules by shocks, as in the case of AGB stars (Cherchneff 2006, 2012), 
and/or the existence of a strong UV field. The great molecular richness of VY 
CMa was already emphasised by Ziurys et al. (2007) with the detection of many 
molecular rotational transitions tracing various components of the wind. Small 
metal oxide species including AlO, TiO, and Ti02 have also been detected in 
the gas-phase close to the star ( [Tenenbaum &: Ziurys 20091 Kamihski et al. 2012, 
2013), revealing a pool of gas-phase species involved in the formation of dust such 
as alumina and titanates. A similarly rich molecular outfiow is found for the RSG 
NML Cygni, but the MOLsphere of a Ori shows a more limited molecular sample, 
including CO, H2O and OH ( |Teyssier et al. 2012D . No SiO lines were detected in 
the Herschel spectra of a Ori reported by Teyssier et al. (2012), in contrast with 
the modelling of spectro-interferometric data by Perrin et al. (2007), for which a 
large SiO column density in the MOLsphere is required. SiO being a key species in 
the formation pathway to silicates, some transitions of this species will probably 
be detected in the wind of a Ori as more Herschel data become available. The 
molecular species so far observed in a few RSGs are summarised in Table 1. 

This plethora of detected molecules indicates that RSG winds must be con- 



Dust Formation in Red Supergiants 



3 



Table 1. Observed molecules in the inner and intermediate wind of the RSG stars a 
Ori, VY CMa, and NML Cyg. Data are taken from Ziurys et al. (2007), Tenenbaum & 
Ziurys (2009), Teyssier et al. (2012), and Kamihski et al. (2013). 
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a Ori 


NML Cyg 


VY CMa 
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OH 
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HCN 
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CS 
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SiS 
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CN 
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PN 
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AlO 
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TiO 
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TiOs 






X 


NaCl 






X 


HCO+ 






X 



ducive to dust formation. However, the dust formation radius derived from data 
modelling in RSGs spans a wide value range. Danchi et al. (1994) derived a large 
inner radius (~ 40 R©) for the dust shell in a Ori from interferometric data, while 
the dust shell for VY CMa was found close to the star (^ 5 R©). Such a location 
is encountered in AGB stars, where dust forms in the shocked gas layers above 
the stellar photosphere ( Cherchneff 2006| . The dust composition in RSGs is still 
a matter of debate. ISO-SWS mid-IR spectra show bands at 9 /im and 18 /im, 
indicative of the presence of amorphous silicates, namely olivine. In the case of a 
Ori, an extra source of mid-IR opacity is required to explain the observed mid-IR 
excess, and the best fit to the data is achieved when alumina (AI2O3) dust is added 
to the composition of the MOLsphere (|Verhoelst et al. 20061 IPerrin et al. 2007)) . 
For a large RSG sample. Speck et al. (2000) and Verhoelst et al. (2006) show 
that the mid-IR excess can be fitted with a combination of 0-rich dust, includ- 
ing olivine, Mgo.8Fei.2Si04, melilite, Ca2Al207, and AI2O3. However, this choice 
of condensates is guided by condensation sequences derived from thermodynamic 
equilibrium (TE) and may then be questionable, as discussed in iJ31 



3 Dust formation as a non-equilibrium process in stellar outflows 

There exists to date no satisfactory model of dust formation in ECEs. The 
great difficulty resides in the characterisation of the different steps involved in 
the formation of solid grains out of the gas phase. In this regard, modelling 
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the production of cosmic carbon dust took advantage of the numerous infor- 
mation provided by the field of combustion chemistry on the formation of soot 
( [Frenklach fc Feigelson 1989[ICherchnefF et al. 19 92). As to silicate dust, notwith- 
standing the studies purveyed by the field of material science on ceramic and silicon 
film growth, the steps leading to the synthesis of silica and silicates are still partly 
unknown. 



In this context and for decades, studies on the synthesis of cosmic dust in vari- 
ous types of ECEs (e.g., AGB stars, RSGs, supernovae) have been based on the as- 
sumption of TE to predict " condensation sequences" (jGrossman fc Larimer 19741 
ITielens 1990l ITielens et al. 1998| ILodders 20061) . and the use of Classical Nucle- 
ation Theory (CNT) to derive grain size distributions for specific dust chemical 
compositions (|Kozasa et al. 19891 |Gail fc Sedlmayr 1999| ITodini fc Ferrara 20011 
IBianchi fc Schneider 20071) . Fo'" 0-rich environments, the condensation sequence 
derived by Tielens (1990) and Tielens et al. (1998) is based on a description of 
the wind assuming TE from the photosphere to large radii, an assumption that 
does not hold in general in ECEs, owing to the complex dynamics of the out- 
flows that include shocks, chromospheric activity, UV fields, etc. For example, 
alumina small clusters form after forsterite clusters in the 0-rich zones of su- 
pernovae ejecta because of the radioactivity-induced chemistry specific to these 
environments dSarangi" fc Cherchneff 2013D . This result is in direct contrast with 
the condensation sequences derived for 0-rich environments (e.g., ITielens 1990p . 
and shows that the use of TE predictions can be misleading in the interpretation 
of complex media such as RSG winds. As to CNT, the description of the dust nu- 
cleation phase assumes steady state, optical properties and concepts (i.e., sticking 
coefficient, surface tension) appropriate for the bulk material. These assumptions 
are clearly not satisfactory to describe the formation of small molecular clusters 
entering the nucleation phase of dust synthesis for the conditions pertaining to 
ECEs. |Donn fc Nuth (1985)| were the first to question the use of CNT for the 
description of dust formation in space decades ago. 



More generally, dust grains form out of the gas phase following a two-step pro- 
cess that involves a nucleation and condensation phase (|Cherchneff fc Dwek 2010l 
ICherchneff 2011( |Sarangi fc Cherchneff 20131 ). The nucleation phase consists of 
the formation of molecules and molecular clusters out of the gas phase and is 
controlled by chemical kinetics and the ability of gas-phase species (i.e, atoms, 
molecules, ions) to react with one another to form larger chemical entities. Chem- 
ical reactions, their rate and occurrence, therefore control the growth of clusters 
from the gas phase. The condensation phase produces dust grains through the 
coalescence, surface growth, and coagulation of large molecular clusters. This 
phase is usually described by various formalisms (moments, bin, or Monte-Carlo 
methods) based on collision theory. Both phases are observed and studied in the 
laboratory with many experimental devices and techniques producing dust such 
as flame aerosol technologies and vaporisation experiments. 
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Fig. 1. The structures of a few clusters involved in the synthesis of forsterite dimers 



taken from Goumans & Bromley (2012) (red: oxygen, yellow: silicon, blue: magnesium, 
grey: hydrogen). Cluster 7 corresponds to the entastite dimer (MgSi03)2- 



4 Nucleation of oxygen-rich dust 



The formation of cosmic dust from a gaseous medium requires high gas densities 
and temperatures and sufficient time for the nucleation of small clusters to pro- 
ceed, followed by their coalescence, coagulation, and growth to form dust grains. 
Therefore, the dense, hot zones above the photosphere of evolved stars such as 
AGBs or RSGs are conducive to the production of dust grains. Because the photo- 
sphere of RSGs is rich in oxygen, silicon, magnesium, and aluminium owing to the 
nucleosynthesis in massive stellar cores, the condensates should include calcium- 
magnesium-iron silicates (members of both the olivine and pyroxene family) and 
metal oxides (e.g., alumina AI2O3 and/or spinel Mg2Al204). Despite the wealth 
of studies dedicated to the formation of silicates and ceramics in the laboratory, 
the identification of molecular precursors to dust formation is still a difficult task. 
For silicates, alumina, or spinel, no gas-phase monomers exist from which solids 
grow by direct coalescence or coagulation. Flame aerosol or shock-induced com- 
bustion experiments have indeed highlighted the important role of heterogeneous 
seed clusters in the formation of these solid compounds (jMcMillin et al. 1996| 
IPratsinis 19981 ICatoire et al. 2003|l . 
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4.1 Silicates 

In an astrophysical context, silicate formation has been investigated in the labora- 
tory by several groups, including Rietmeijer et al. (1999, 2002), K"aito et al. (2003)[ 



Kamitsuji et al. (2005) and Kimura et al. (2008) In all experiments, the va- 



porisation of solid magnesium and silicon oxides led to the formation of crys- 
talline Mg-bearing silicate like forsterite, Mg2Si04, Mg-bcaring amorphous sili- 
cates, MgO grains (periclase and amorphous), Si02 grains and pure Si grains. 
Several studies emphasise a two-step process in which periclase and silica form 
in the first place, and further coalesce to give forsterite or amorphous silicates. 



e.g., Kimura et al. (2008) Under circumstellar conditions, |Kohler et al. (1997) 



showed that MgO grains were not the preferred intermediates in the dust nucle- 
ation process, a point corroborated by recent chemical models of supernova ejecta 
( [Sarangi fc Cherchneff 2013 1. In the specific context of H-rich ECEs (e.g., AGB 



and RSG winds), the large H2 content ensures that reaction of II2 with MgO, a 
highly exothermic process, will destroy the newly formed small oxide, hampering 
its formation in the gas (jPlane 2013P . 



Recently, Goumans & Bromley (2012) proposed a new nucleation scheme lead- 



ing to small silicate clusters with the stochiometric composition of enstatite and 
forsterite. The scheme is based on the more energetically favourable paths to sil- 
icate growth. The growth mechanism consists of the initial formation of the SiO 
dimer, Si202, from a SiO-rich gas, followed by a two-step sequence involving 1) 
the addition of one oxygen atom to SiO dimers and subsequent silicate intermedi- 
ates through reactions with O-bearing species available in the gas phase; and 2) 
the addition of one magnesium atom to the silicate cluster lattice. All chemical 
reactions characterising these addition steps are down-hill processes (exothermic) , 
unlike the formation of the SiO dimer. Therefore, the closure of the Si202 ring 
represents the bottleneck to the synthesis of small silicate clusters. Some inter- 
mediate clusters involved in the formation of enstatite and forsterite dimers are 
illustrated in Figure 1. 



4.2 Alumina 

The structure of small alumina clusters has been extensively studied theoreti- 
cally (e.g., lArchibong fc St-Amant 1999[ |Li fc Cheng"20T2l ) , and some structures 
of small AI2O3 polymer clusters are illustrated in Figure 2. Although there exists 
no well-identified chemical pathways to their formation, a possible route can be 
conjectured based on cluster structures. Because of the structural similarities be- 
tween (SiO)2 and (A10)2, the formation of small alumina clusters may first involve 
the closure of the (A10)2 ring, and the later addition of oxygen to lead to AI2O3. 
Further growth to AI2O3 dimers and trimers may occur with further addition 
of Al via reaction with Al, AlO, and AlOH, while oxygen addition may proceed 
from reaction with II2O, OH, and/or O2. Such a scheme awaits confirmation by a 
thourough assessment of the chemical pathways and rates involved in each growth 
step. 
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Fig. 2. Examples of the most stable structmes of small AI2O3 clusters adapted from Li 
& Cheng (2012) (red: oxygen, grey: aluminium). 



5 Silicate clusters in 0-rich AGB winds 

The formation of silicate dust in RSGs chemically resembles that occurring in 
0-rich AGB stars. Both environments are H-rich and comprise silicon monox- 
ide, SiO, atomic magnesium and iron, and water. To test the silicate chemical 
growth proposed by Goumans fc Bromley (2012)[ t he inner wind of the Mira star 
IK Tau is modelled following IDuari et al. (1999)1 and Cherchneff (2006). The 



model includes a chemical kinetic description of the various gas-phase processes 
coupled to the nucleation routes to forsterite dimers, and considers the impact of 
shocks on the outer atmosphere of IK Tau. This approach was successful in repro- 
ducing the chemical composition of the inner wind of carbon- and 0-rich AGBs 
(|Cherchneff 2006llCherchnefF 2012"!) . The abundances of small silicate clusters with 
respect to II2 as a function of pulsation phase in the post-shock gas above the pho- 
tosphere are illustrated in Figure 3. One striking aspect is the rapid conversion of 
SiO to SiO dimers and the growth of clusters leading to the formation of forsterite 
and cnstatite dimers when the post-shock gas temperature drops below ~ 2500 
K. This formation does not equally deplete all clusters, and small clusters such as 
SiO dimers and Mg3Si207 may survive the nucleation step. These clusters may 
either be expelled as free-flying species in the wind or eventually coalesce at lower 
temperatures to form amorphous silicates. This preliminary study indicates that 
silicate dust clusters seem to easily form in the inner wind of Mira stars, pointing 
to the efficient formation of silicate dust in these environments. 



6 Conclusions 



The gas conditions in RSG winds above the photosphere are hasher than those 
found in AGB outflows. As revealed by optical and UV observations of a Ori, the 
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0.4 o.e 

Phase 



Fig. 3. Abundances of silicate clusters with respect to H2 as a function of pulsation 
phase in the shocked inner wind of the Mira star IK Tan (Gobrecht & Cherchneff, in 
preparation) . 



chromospheric activity generates a hot plasma that coexists with a more abundant, 
cooler gas lifted above the photosphere by convective cells (|Lim et al. 1998p . It is 
then likely that the formation of dust in RSGs occurs in the dense, cool, clumpy 
component of the wind. A puzzling case is that of a Ori, where AI2O3 dust is 
required in the MOLsphere close to the star to account for the mid-IR opacity 
excess, whilst silicate grains are observed at large radii above the stellar surface. 
Two scenarios have been proposed. The first scenario considers a two-shell con- 
figuration where dust destruction by chrosmospheric activity is effective after the 
formation of AI2O3 until large radii and cooler temperatures are reached where 
anew silicate dust forms (jVerhoelst et al. 2009[) . A second scenario considers the 
steady formation of AI2O3 close to the photosphere. AI2O3 being transparent, it 
would become detectable at radii where silicates form and coat the AI2O3 grains, 
resulting in a silicate dust shell at large radii (|Onaka et al. 1989)) . Both scenar- 
ios overlook the problem of the low gas density encountered at large radii which 
will hamper efficient silicate cluster coalescence and coagulation at these large 
distances. 

In some RSGs like VY CMa, the dust formation process may resemble that 
operating in AGB winds, where molecules and dust clusters are efficiently produced 



Dust Formation in Red Supergiants 



9 



close to the star. For these objects, dust may contribute to the triggering of the 
mass loss close to the stellar surface. In the case of a Ori, the low molecular 
content and poor chemical complexity observed close to the stellar photosphere 
(see Table 1) may indicate a less efficient dust production process. Dust there 
may not be a key player in the wind acceleration mechanism. More generally, dust 
formation is not an equilibrium process and is highly dependent of the physics and 
chemistry of the locus where grains form. Because the dust formation zone in RSGs 
is clumpy and has a complex dynamics, it is difficult to derive plausible scenarios 
for the formation of dust in such environments without a proper physico-chemical 
model of RSG winds. 
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